Effects of coupling to vibrational modes on the ac conductance of molecular junctions 
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We theoretically examine the effect of the coupling of the transport electrons to a vibrational mode of the 
molecule on the ac linear-response conductance of molecular junctions. Representing the molecule by a single 
electronic state, we find that at very low temperatures the frequency-dependent conductance is mainly enhanced 
(suppressed) by the electron-vibration interaction when the chemical potential is below (above) the energy of 
that state. The vertex corrections of the electron-vibration interaction induce an additional peak structure in the 
conductance, which can be observed by tuning the tunnel couplings with the leads. 

PACS numbers: 71.38.-k, 73.63.Kv, 73.21. La 



I. INTRODUCTION 

The fabrication of junctions made of small molecules seems 
by now well established. 112 The transport properties of such 
molecular junctions are largely determined by the interplay 
between the electrical and the vibrational degrees of freedom. 
These include the couplings of the bridge with the leads, as 
well as the resonance energies. Their effects on the dc conduc- 
tance of such junctions have been examined in many works 
(see below). Here we study the ac conductance, concentrating 
on the modulation of the frequency dependence brought about 
by the coupling to the vibrational modes. Our calculation is 
confined to the linear-response regime, and is carried out to 
second order in the coupling with the vibrational modes. 

The simplest model for describing this system assigns a sin- 
gle electronic level, e , to the molecule, which is connected 
via two leads to two electronic reservoirs. Those are kept at 
different chemical potentials. When electrons pass through 
the molecule, they are coupled to its vibrational modes. Since 
the molecule has a finite size, those modes have finite frequen- 
cies; in the simplest approach and at very low temperatures, 
they may be represented by a single vibrational frequency, 
w , treated within the Einstein model. Various approximate 
numerical and analytical schemes have been developed for 
treating this model. The dc current-voltage characteristics 
have been studied by the perturbation theory^— the rate 
equations^ and Monte Carlo methods. 21 Furthermore, the 
shot noise and the full counting statistics of the charge pass- 
ing the junction have been studied by perturbation theory* 2 ^— 



i-H , rate equations, and the polaron approximation 
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In the linear-response regime, the Breit-Wigner resonance 
of the dc conductance (as a function of the equilibrium chem- 
ical potential fi, and at very low temperatures) is narrowed 
down by the electron-vibration (e-v) interaction, due to the 
renormalization of the tunnel coupling between the molecule 
and the leads (the Frank-Condon blockade) However, the 
e-v interaction does not yield side peaks in the conductance at 
energies corresponding to the molecule vibrational frequen- 
cies, as long as it is not large enough to induce polaronic 
phase transitions. For such peaks to arise (at zero temperature 
and at small values of the e-v coupling), the electron needs 
to lose an energy w (we use ft = 1) in order to excite a vi- 
brational mode, and this is not possible in the linear-response 
regime. Such side peaks will appear (at zero temperature) in 



the nonlinear-response regime, where the finite voltage, V, 
allows for the opening of inelastic channels when eV exceeds 
uj JdiI2. Side peaks can be induced by a finite temperature, but 
their height will be minute, reflecting the number of available 
vibrational modes.— 

The ac conductance of molecular junctions, in particu- 
lar the effects of the e-v interaction on its frequency depen- 
dence, has been studied to lesser extent. Recently, Kubala and 
Marquardt 30 wrote down expressions for the ac conductance 
of interacting electrons traversing molecular junctions, in the 
case where the dc voltage is finite, while the ac one is very 
small, including as an example the e-v interaction. Our ex- 
pressions (which are derived by a different approach) agree 
with theirs in the zero-dc bias limit. In addition, we give de- 
tails of an extra diagram [(d) in Fig. Q~| which they chose to 
ignore (because it vanishes in the dc limit), and present a de- 
tailed discussion of the dependence of the ac conductance on 
the various parameters. 

As we show, the e-v interaction generates four contribu- 
tions to the ac conductance. The first two are a Hartree and an 
exchange terms. In addition, the Kubo formula for the ac con- 
ductance also contains two vertex corrections. In the Keldysh 
formalism, 31 - 32 the vertex corrections correspond to an expan- 
sion of the self energies in the time-dependent chemical po- 
tentials of the left and right reservoirs, S/j, l (t) and Sfi R (t), re- 
spectively. These corrections cause the ac conductance to di- 
verge when the frequency u of the applied ac voltage crosses 
lu q ; the imaginary part of the Green function of the vibrational 
mode is infinitesimally small to second-order in the e-v cou- 
pling. To avoid this divergence, a finite lifetime of the vibra- 
tional modes should be included. Two effects can be consid- 
ered which contribute to this lifetime. One is the relaxation 
due to the coupling with bulk phonons in the substrate. We 
denote this relaxation rate by S; the other involves the elec- 
trical polarization, which results from the coupling of the vi- 
brational mode with the transport electrons. The latter effect 
is dominant for floating molecules, air-bridged between the 
electrodes (in the absence of the substrate)^ We discuss this 
case in Sec. |lV]and in Appendix 151 

The contributions of the vertex corrections to the ac conduc- 
tance are proportional to the combination T L 5fi L + r R 6/i R , 
where r L ( R ^ denotes the broadening of the energy level 
representing the molecule due to its coupling with the left 
(right) lead. Choosing a symmetric configuration, Sfi^t) — 
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—6fi K (t) and T L — T R , causes these contributions to van- 
ish. Therefore, in this symmetric case the ac conductance is 
affected by the e-v interaction only via the Hartree and the 
exchange (Fock) contributions. As we show, the former con- 
tribution is usually dominant, causing a relatively small in- 
crease (larger decrease) of the ac conductance when the av- 
erage chemical potential in the leads, \i, is below (above) the 
electronic level e on the molecule. Similar to the dc conduc- 
tance, this causes a narrowing of the ac conductance (plotted 
versus /i). Since the exchange contribution is even in (/i — e ), 
subtracting the conductances above and below e can yield in- 
formation on the Hartree contribution. When the junction is 
not fully symmetric, the vertex corrections introduce an addi- 
tional structure in the frequency dependence of the ac conduc- 
tance, especially near ui = ±w . These corrections become 



maximal in the limit I\ = V and T Y 



0. 



The organization of the paper is as follows. We begin in 
Sec. |II]by describing our model and presenting the expression 
for the ac conductance of the system. The details of the deriva- 
tion are relegated to Appendix|A] SectionQJJJis devoted to the 
analysis of the results. In Sec. [IV] we briefly discuss pro- 
posals for possible measurements and the effect of the elec- 
tronic polarization on the conductance. The random-phase- 
approximation (RPA) treatment of the latter effect is outlined 
in Appendix IE1 



II. THE ac CONDUCTANCE 

Our model system consists of two electronic reservoirs, 
connected together via a single electronic level e , which rep- 
resents the molecule. The left and right reservoirs are kept 
at time-dependent chemical potentials, /x L (i) = fj, + Sfi h (t) 
and n R (t) = fi + <5/i R (t), which oscillate with frequency lu. 
When the electron is on the molecule, it is coupled to a local 
vibrational mode of frequency uj . This simplified model is 
described by the Hamiltonian 



(1) 
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It consists of the leads' Hamiltonian 

k p 

the Hamiltonian of the molecule 

U mo \ = £ c o c o + u Q tfb + 7(6 + tf)c\c Q , (3) 

and the tunneling Hamiltonian describing the coupling be- 
tween the molecule and the leads 

Htu„ = $>l4co + h.c.) + £>r4 c o + h.c). (4) 

k p 

Here, cL s and c fc , p % denote the creation and annihilation op- 
erators of an electron of momentum k(p) in the left (right) 
lead, respectively. The creation and annihilation operators on 
the level e are denoted by cj and c , and 6^(6) creates (an- 
nihilates) a vibrational mode of frequency u> . The coupling 



of the transport electrons with the vibrational mode is scaled 
by 7. The broadening of the resonant level on the molecule, 
T = T L +T R , is given by r L (R) = 27w|£ L( . R ^ | 2 , with v being 
the density of states of the electrons in the leads. 

The current flowing into the molecule from the left reser- 
voir may be expressed in terms of the Keldysh Green 
functions, 

J L (t) = -2eRc f dt>Y,\tL\ 2 [G r 00 (t,t')g<(t>,t) 
J k 

+ G<(t, t')g a k (t' -t)- g r k (t - t')G<(t>, t) 
-g<(t,t')G a Q0 (t',t)}, (5) 

with an analogous expression for the current coming from the 
right reservoir (with L replaced by R and k replaced by p). 
Here, 

G r 00 (t,t') = - i 9(t-t / )({c (t),cl(t')}), 
G< (t,t')= i (c (t / )cl(t)}, (6) 

are the retarded and lesser Green functions on the molecule, 
where (• • • ) denotes a quantum and statistical average over 
the states of the whole system. On the other hand, the retarded 
and lesser Green functions on the left (right) lead, 



g r k{p) (t,t') = -iO{t-t')({c k(p) {t) lC l {p) {t')}) a , 



(7) 



are given by the average (• • • ) over the noninteracting leads 
alone. 

FigureQ]depicts the diagrams of the linear-response ac con- 
ductance: (a) in the absence of the e-v interaction, (b) includ- 
ing the self energy from the Hartree term and (c) including 
the exchange term. Diagrams (d) and (e) are the vertex cor- 
rections of (b) and (c), respectively. The solid line indicates 
the Green function of the electron on the molecule whereas 
the dotted line denotes the Green function of the vibrational 
mode. The wavy line indicates the frequency of the ac field. 
The diagrams in Fig. Q]are a ls° listed in Ref. |3(J However, 
the diagram (d) was not calculated in that reference. Our ex- 
pressions for the other diagrams agree with those of Ref. 13CJ 

We expand Eq. (O to first order in Sfi L (t) and (5/i R (t) as de- 
tailed in Appendix lAl The linear-order currents if go ( a) ) are 
obtained by substituting the expansions ( lAlOt and d A 1 1 b into 
Eqs. ( 1A1I ) and JA21 I. and then substituting the latter expan- 
sions into Eq. (O. In our two-lead system, the linear-response 
ac conductance can be expressed as a matrix, 



(8) 



We show in Sec. [Till that when the molecule is coupled sym- 
metrically to the leads, T L = r R , and the ac voltage is applied 
symmetrically as well, <5/i L (w) = — (5/z R (w) = 5n(ui)/2, the 
conductance pertaining to the net current is given by 
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FIG. 1: The diagrams of the ac conductance: (a) without the e- 
v interaction, (b) the Hartree term, (c) the exchange term, (d) the 
vertex correction for the Hartree term and (e) the vertex correction 
for the exchange term. The solid line denotes the Green function of 
the electrons while the dotted line indicates the Green function of 
the vibrational mode. The wavy line indicates the frequency of the 
external ac field. 



In this case, diagrams (d) and (e) (see Fig. [TJ vanish since 
A 2L = A 2R = [seeEqs. dA16); A 2L(r) = r L(R) (r L <5/x L + 
r R <5/i R )]. Consequently, the conductance given in Eq. ©can 
be presented as a sum of three terms, 



= G nint M + Gh(w) + G«(w), 



where 



G „intM = e 

G H (w) = e 



J R-ni„tM 



25n(u) 
26 

7 L-ex( W ) - 7 R-exM 



2S/x(w) 



and the partial currents appearing in Eqs. (fTTl i are given in 
AppendixlAl see Eqs. ( IA15I ), dA19t . and (IA21b . The conduc- 
tance pertaining to the fully-symmetric junction, [see Eqs. (O 
and (TTOb l is analyzed in the next section. 

In the maximally asymmetric tunneling configuration, 
where T L = T and T R = 0, the conductance of the junction 
is given by 



(12) 



In this case there are contributions also from the vertex cor- 
rections, diagrams (d) and (e) of Fig. Q] We analyze the con- 
ductance of this configuration in Sec. [Ell] expressing it in the 
form 

G(u>) =G nint M + G H M 

+ G 0X H + G vcrH H + G (w). (13) 

Here, each of the partial conductances is given by the respec- 
tive partial current [see Eqs. (|A15t , dA19t , ( IA2U . (IA221 >. and 
(lA23l l, divided by Sfi h (u))/e. 
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FIG. 2: The ac conductance of a fully-symmetric bridge, as a func- 
tion of e — fi at u> = r. The coupling strength of the e-v interaction 
is 7 = 0.3T, and the frequency of the vibrational mode is ui Q — T 
(dashed line) and uj = 0.5T (dotted line). The solid line is the 'bare' 



conductance G„ 



obtained in the absence of the e-v interaction. 



HI. NUMERICAL RESULTS FOR THE ac CONDUCTANCE 

As explained in Sec. U one has to allow for a finite lifetime 
(10) for the vibrational mode (even at lowest-order in the coupling 
7), in order to avoid unphysical divergences. In this section 
we assume a constant relaxation rate (adopting the value 5 = 
O.ir for the numerical computations); possible contributions 
of the transport electrons to this rate are discussed in the next 
section. In the following we measure all energies in units of 
r, the broadening of the resonance molecular electronic level 
brought about by the coupling to the leads. 

We begin with the case of a symmetrically-coupled junc- 
tion, T L = r R . Figure [2] shows the conductance G, Eq. (O, 
(ID as a function of e — fx, for a fixed ac frequency uj = T. The 
solid line indicates G n!nt , the conductance in the absence of 



T nint 

the e-v coupling. It is seen that the resonance peak does not 
reach unity as happens in the case of the dc conductancer^^ 
due to the suppression by the ac field. When the e-v coupling 
is accounted for (the dashed and dotted lines in Fig. |2]i this 
peak becomes somewhat narrower, which is more remarkable 
for smaller lj q . This behavior of the resonance peak is quali- 
tatively similar to the dc case^S 

Furthermore, the center of the resonance peak shifts to a 
higher energy as compared to the one in the absence of the 
e-v coupling. The narrowing and the shift may be understood 
by including the Hartree term self-consistently in the electron 
Green function [Eq. ( 1A12H employing the Dyson equation 



ea/[ w - Eo + »r/2-£&(o)]. 



(14) 



As can be seen from Eq. ( 1A20I ). S R (0) is negative, and conse- 
quently the localized level representing the molecule is shifted 
to a lower energy, implying a lower (higher) conductance for 
e < /i (e > /i). Also the Hartree term renormalizes the 
energy scale, which contributes to the narrowing. This renor- 
malization is due to the polaron binding energy. Notice that 
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FIG. 3: (a) The ac conductance of the fully-symmetric bridge, as 
a function of the frequency of the bias voltage cj, at Sq — /J, = T 
and 7 = 0.3r. The frequency of the vibrational mode is cu = T 
(dashed line) and uu = 0.5r (dotted line). The solid line is the 'bare' 
conductance G nint , in the absence of the e-v interaction, (b) The 
additional conductance due to the Hartree term of the e-v interaction, 
(c) The additional conductance due to the exchange term. 



the energy scale is also renormalized by the exchange term, 
which oppositely tends to broaden the resonance peak. How- 
ever, this effect is smaller than that of the Hartree one (as long 
as 8 remains small). Since in the linear-response regime there 
is no real exchange of energy between the electrons and the 
vibrational mode, no additional peak structure appears in the 
conductance. 

The full conductance G, Eq. (O, as a function of the ac 
frequency u of the external bias when e Q — fj, = V, is depicted 
in Fig. [3] (a). The solid line indicates the conductance in the 
absence of the e-v coupling, G nint . Two broad peaks appear 
around u of order ±1.5(£ — fi). The broken lines show G in 
the presence of the e-v interaction with oj a = V or lu q = 0.5F. 
The e-v interaction increases the conductance in the region 
between the original peaks, shifting these peaks to lower \u\, 
while decreasing the conductance slightly outside this region, 
where the e-v effect decays very quickly. Figures [3] (b) and 
(c) portray the contributions to the conductance due to the e-v 
interaction coming from the Hartree and the exchange terms 
[diagrams (b) and (c) in Fig. Q~|, G H and G cx respectively, 
for the same parameters. Similar results arise for all posi- 
tive e — /i. Both G H and G cx show two sharp peaks around 
uj ~ ±(e — fi) (causing the increase in G and the shift in 



its peaks), and both decay rather fast outside this region. In 
addition, G ex also exhibits two negative minima, leading to 
small 'shoulders' in the total G which are not much visible in 
Fig. [3(a). The exchange term virtually shifts the polaron level 
on the molecule, yielding an enhancement in the conductance. 
The amount of increase is more dominant for lower u> . The 
situation reverses for e Q < /i, as seen in Fig. [4] Here, G nint 
remains as before, but the ac conductance is suppressed by 
the e-v interaction. The additional conductance is dominated 
by G H . The Hartree term of the e-v interaction renormalizes 
the energy level in the molecule to lower values, resulting in 
the suppression of G, while G cx coincides with the case for 
e Q — [a = T. The amount of decrease is larger for lower ui . 
Figures 3 and 4 show that the sign of change in the conduc- 
tance depends on the sign of the Hartree contribution, Gh • 

Next, we consider the conductance Eq. ( fT3l ) of an 
asymmetrically-coupled bridge for which T L = T and T R = 
0. In this case the vertex corrections are required as explained 
in Sec. [TT] Figure|5]shows the conductance as a function of the 
bias frequency u> when e — /.i = T (panel a) and e — fj, = — T 
(panel b). Similar to the conductance of a symmetrically- 
coupled junction (T L = T R ), the ac conductance is enhanced 
or suppressed as compared to the noninteracting case, depend- 
ing on whether e > \i or e < /i. The contributions to the 
conductance arising from the vertex corrections to the Hartree 
and the exchange diagrams, G vcr H and G VG rex [see Eq. ( fT3] ll 
are shown in panels (c) and (d) respectively. The vertex cor- 
rections coincide for e — fj, = ±T. Interestingly, the plot (c) 
exhibits sharp peaks at lu — ±w due to the singularities of the 
Green function of the vibrational mode, which are smeared by 
its lifetime. These also appear as anomalous peaks in the to- 
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FIG. 5: The ac conductance of an asymmetric junction, as a function 
of oj at (a) e — fi = T and (b) e — n — — T when T L = T. The 
energy of the vibrational mode is ui Q = T (dashed line) or tj = 0.5r 
(dotted line). The solid line is the 'bare' conductance G nint , in the 
absence of the e-v interaction. The additional conductances resulting 
from the vertex corrections of the Hartree and exchange terms are 
shown in panels (c) and (d). 



tal conductance. This is due to the charge fluctuation caused 
by the ac field, which leads to a fluctuation of the vibrational 
mode. This effect becomes blurred when the relaxation rate 5 
is large. Note that the contributions of the vertex correction of 
the exchange diagram, G vorox , cancels some amount of G cx , 
therefore the enhancement of the conductance is not large for 
e > fj,. 



IV. DISCUSSION AND SUMMARY 

Figures[3]and|4]show that within our perturbative expansion 
the effect of the e-v interaction on the ac conductance is rela- 
tively small and the results do not differ qualitatively from the 
expectation without this interaction. One way to focus on the 
e-v effect is to compare the conductances at e — fj, = T and 
e — /i = — r, As mentioned in Sec. [TTTJ all the contributions 
to the conductance except the one coming from the Hartree 
term, G H , are the same for e — /i = T and e — \i = —F. 
Therefore, one can measure the difference G H — G' H by sub- 
tracting the total conductance G' for e — /i = —Y from G 
for e — /i = r. Another effect which can be observable 
involves the sharp peaks due to the vertex correction of the 
Hartree term. This effect is only visible when the combina- 
tion r L (5/i L +r R (5// R is large and the relaxation time 5 is small 
compared to T. Tuning the tunneling energies r L ( R ) can help 
to extract these effects from the data. 

As stated, the vibrational lifetime for floating molecules 
(that are not placed on a substrate) may be dominated by the 
effect of the electronic polarization (i.e., the self energy of 
the vibration due to the coupling with the transport electrons), 
which also smears the singularities of the Green functions of 
the vibrational mode D a (a takes the values r, a, or <) in Eqs. 
( |A6t and ( 1A7| >. In this case, at order r ) 2 , D a is replaced by 
Eqs. (ED 

or (IB2I > in Appendix|B] which include the RPA-type 
dressing by electrons shown as the thick dotted line in Fig. [6] 
Thus, all the dotted lines in Fig. Q]are replaced by the thicker 
dotted line. As discussed in Appendix 151 this procedure gen- 
erates an additional lifetime for the vibrational modes. For 
7 2 /(w r) < 1/4 we find that the effect of these RPA-type 
dressings on the conductance is similar to that of the effective 
'ad hoc' relaxation rate 6, as reported previously. 33 For larger 
e-v coupling 7, the energy shift of the electronic level repre- 
senting the molecule e [which corresponds to E^(0) in Eq. 
(IA18H is no longer monotonic as a function of e . This may 
imply the breakdown of the perturbation expansion. 

As mentioned in the introduction, and seen in Eqs. ( IA22I ) 
and dA231 >. the vertex contributions [diagrams (d) and (e) of 
Fig. Q~) appear only for asymmetrically-coupled bridges. It is 
diagram (d) of Fig. Q]that is most sensitive to the vibrational 
mode lifetime, since its contribution is explicitly proportional 
to the vibrational Green function (the dotted line in the cen- 
ter of the diagram). Our results for the symmetric case, based 
on diagrams (a)-(c), are indeed not very sensitive to the vibra- 
tional mode lifetime. The RPA treatment of the electronic po- 
larization for the case in which the molecule is coupled asym- 
metrically to the leads requires several additional terms. In 
addition to diagrams (d) and (e) of Fig. [T] in the asymmet- 
ric case we also need to add diagrams (f) and (g) of Fig. [6] 
which represent the vertex corrections to diagrams (b) and (c) 
of Fig. Q] However, these new contributions involve higher 
orders of the e-v coupling 7, and thus go beyond the scope of 
the present paper. 

In summary, we have studied the effect of the e-v inter- 
action on the ac linear conductance of molecular junctions, 
employing a simple model. The e-v interaction enhances or 
suppresses the conductance, depending on whether the energy 
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level of the orbital in the molecule is higher or lower than the 
chemical potential. When the tunnel coupling is asymmetric 
as defined in Sec. [Il]above, an additional anomalous structure 
appears at u) = ±w due to fluctuations caused by the ac field. 
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Appendix A: Details of the current derivation 

As mentioned, our calculation is carried out to second-order 
in the coupling of the transport electrons to the vibrations, 
7. Thus the Dyson equations for the Green functions on the 
molecule are^i 



G r ^{t,t') = G r ^\t-t') 

+ J d^G^^Ci-t^S^^^JG^^^-O + J dt x J ^ 2 G^ K0) (t-t 1 )S e t ] (t 1 ,* 2 )G^ ](0) (t 2 -O ) (AD 



G< (t,t') = G<°\t,t') + / dt 1 G^\t-t 1 )Z r H (t 1 ,t 1 )G< i0) (t 1 ,t') + J dt 1 G^°\t,t 1 )^ l \(t 1 ,t 1 )G^°\t 1 -t') 
+ J dt x J dt 2 G^ ) (t-t 1 )EUt 1 ,t 2 )G< 0) (t 2 ,t') + J dt x J dt 2 G^\t,t 1 )^ K {t 1 ,t 2 )G a ^ ) {t 2 -t') 
+ J dt x J dt 2 G r f{t - - t'). (A2) 



Here, T,\\(t,t) represents the self energy due to the Hartree 
term [diagram (b) in Fig. [T|, 

E£(t,t) = -i! 2 I dt'G< {Q \t',t')D r (t - 0, (A3) 



and t) — £jj(t,t). The contribution of the Fock term 

[diagram (c) in Fig. Q~|, £g X (i, 0> wnere a takes the values r, 
a or <, is 

Zi(t,t') =i 1 2 [G<°\t,t')D a (t-t') 

+ G&°\t - t')D<(t - ± G^°\t - t')D a (t - t% 

(A4) 



and 



£< (*,*') = i 1 2 G< {0) (t,t')D < (t - t'). (A5) 



Here, G^ 0) (t-t') are the Green functions in the absence of 
the e-v interaction. 

The Green functions pertaining to the vibrational mode are 
given by 

D r (t - = -iO(t - t')([b(t) + &t(t)], [&(*') + ft* (*')]) 

(A6) 



and 



D<(t - = -i([b(t') + b\t% [b(t) + 6t(t)]>. (A7) 

To second-order in 7, these are required only to order 7 . 
However, as explained in Sec. |U there is a need to assign a 
finite lifetime to the vibrations. Here we take the relaxation 
rate (i.e., the inverse lifetime) to be a constant, 5, arising from 
a possible coupling to a substrate. A different scenario is con- 
sidered in Appendix [B] Thus, the Fourier transforms of the 
vibrational Green functions (at zero temperature) are 



and 



1 



1 



uj — uir, + iS uj + uj n + iS 



(A8) 



D<H = 



uj + u! n + iS uj + ui n 



iS 



(A9) 



Note that all the Green functions except Gqq ^ (t, t') and 

ffoo (*! depend only on the time difference t — t' . 

We expand the Green functions pertaining to the elec- 
trons to the first order in the chemical potentials S/j, l (uj) and 



7 



9< {p) (t,t') « i/[& (p )]e-«MP)(*-*') exp[i J* dhS^it)] = i/fo w ]e-'W t - t '> {1 - | ^[e~^ - *> Ml(r) (c)}, 

(A10) 

and 

(&®M<*f ^e-*-(*-o i/(w) rGS(°)( w )GSS 0, («) 

+ ^'^^[/(^W)-/^)^^^ (All) 



where £, = £, — and /(2) = l/[exp(w/fcBT) + 1] is the 
Fermi distribution function of the leads. The Fourier trans- 
form of the retarded Green function of the molecule is 



1 



072' 



(A12) 



and the advanced Green function is G 



a(0) 



M = [c^ 0) Mr 



Once the current emerging from the left lead, Eq. ©, is 
expanded to linear order in S/j, l and <5/i R , it is convenient to 



divide its Fourier transform into five parts, 

+ ^L-vcr H H + ^L-vcrcxM- (A 13) 



The first term on the right-hand side of Eq. dA13l ) is the cur- 
rent in the absence of the e-v interaction, 



J L-nintM 



2lTUJ 



du>'\f(w + w') - /(a;0]Re{-iA 1L [G^ O) (a; + «') - G&V)] + A 2L G^ 0) ( W + u/)G& 0) (a/)} . 

(A14) 



Carrying out the frequency integration, we find 

^_ nint (2) = ^Re([(£5 - 1)A 1L 



1 



22(i + cj_ 

+ log[l + 4(2 - e ) 2 ] + log[l + 4(2 + e ) 2 ] - 2 log[l + 4eg)]}) 

I 



{— 2i arctan[2(cj — e )] — 2i arctan[2(cj + e ) 



(A15) 



We have introduced in Eq. (IA14b the notations 

A 2L = r L( r L% + r R^R)- 



(A16) 



In Eq. JA15t and below, we denote by tilde all variables di- 
vided by r (assigning for them the same notations as before). 

The second term on the right-hand side of Eq. JA13I ) is the 
contribution of the self energy resulting from the Hartree term 
E H [diagram (b) in Fig. [H, 

Jl-hM = T~ I du'[f(u + u') - /(«')] 



+ A 2L {[c4VW)] 2 GjV) 
+ G^ ) ( W + c')[G^ 0) (^)] 2 })^h( 



,^12^0(0), 



(A17) 



Here, EJj(O) is 
S&(0) 



r 

2?r 



duT(?P\ U )G$ s >(cj)D r (0). (A18) 



Performing the integration in Eq. (IA17b yields 



^_ H (2) = ^Re 



32e [(i2 - 1)A 



1L 



1 2LJ 



(l + 4g2)[4^-(i + 22)2 



with 



E£(0) = - — [tt - 2arctan(2e )]: 



(5 2 



e&(o), 

(A19) 



(A20) 
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The contribution of the exchange terms [diagram (c) in Fig. [TJ to the current is 



2ituj 



dw'[/(a, + u/) - /( W ')]Re(-zA 1L [G^ 0) ( W + u/)££> + u')G r ™(u> + 



- G o 1%0^')G n 1°V)] + A 2L {[G^% + w')] 2 K x (" + "')<j°V) + <4% + c/)[G^(w')] 2 5- (u/)}). 

(A21) 

Finally, the last two terms on the right-hand side of Eq. (1A131 > are 



and 



^L-vcrH^) — 



e7 



(2tt) 2 w 



du/ / rte{ J A 2L [/( W | W ')-/( W ')]G;f(, + ,')GS U) ( W yH 



/^a(0)/ n 



Re 



i 2L 



(2tt) 2 \25(i + 5) ; 



r{— 2i arctan[2(o; — eq)] 



2i arctan[2(5 + ?„)] + log[l + 4(5 - e ) 2 ] + log[l + 4(5 + ? ) 2 ] - 2 log[l + 4£§]}'- g° - 

(w + io) 2 - 



e7 



(2tt) 2 w 



dw' / dw"Re{zA 2L [-/(w + u/)L> a (u/ - w") + f(uj')D r (uj' - u") 



+ D<(u,> - + w") - /(w")]G^ 0) (w + ^)G tV)G[£V + w")GSo°V)}- 

I 



(A22) 



(A23) 



i L-ex 



We evaluate the integrals appearing in the expressions for 
-^L-verex [Eqs. (IA2U and (IA23H numerically as 
functions of w, u Q , and e . The linear expansion of the current 
emerging from the right lead, ijj, is calculated in the same 
way. 



Appendix B: Green function of the vibrational mode including 
self energy corrections due to the transport electrons 

Here we solve the Green function of the vibrational mode 
taking into account self energy corrections coming from the 
coupling with the transport electrons. We derive the elec- 
tronic polarization induced by those electrons employing the 
random-phase-approximation (RPA). 

The relevant diagrams are depicted in Fig. [6] leading to 

dtj dt 2 D^°\t-t 1 )W^{t 1 ,t 2 )D^ a \t 2 X), 

(Bl) 

and 

D<(t,t') = Jdt.J dt a ir(t,t 1 )n<(t 1 ,i a )D o (t a ,0. 

(B2) 

Here, D a [see Eqs. <|A6) and (lATt l and D a( ~^ are the Green 
function pertaining to the vibrational mode in the presence 
and absence of the coupling to the transport electrons, respec- 
tively; n Q is the polarization, 

+ G^(t-t')G<°\t',t)], (B3) 



and 



n<(t,t') = -i^G<°\t,t')G>^(t',t). 



(B4) 



The expressions for the polarization H a contain G <i ~ \ which 
depends on the chemical potentials in the reservoirs. Expand- 
ing to linear order in <5/ij/ R s, yields diagrams (f) and (g) in 
Fig. [6] This order of the expansion turns out to be proportional 
to T L (5/i L + T R (5/z R , and therefore diagrams (f) and (g) of Fig. 
[6] do not contribute to the conductance of a fully-symmetric 
junction [Eq. ©J. In any event, these diagrams involve higher 
orders in 7, and therefore we discard them hereafter. 

Confining ourselves for simplicity to the case of a fully- 
symmetric junction [see Eq. ©], it is sufficient to keep only 
the zeroth order for G < (°\ which yields 



1 



2uj q 



and 



where 



r5 2 -5 2 - 2S ir-M(w)/r 



D < (u) = L> r (5)n<(5)L> a (5), 



(B5) 



(B6) 



IT(0) 



27 2 r 



when 5 = 0. When w^Owe have 



(B7) 
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7 2 r r i 

n r (w) = —i — {2i arctan[2(w — e )J + 2i arctan[2(w + e ) 

2tt 2i + 2u> 

+ 2 ln[l + 4eg] - ln[l + 4(w - e ) 2 ] - ln[l + 4(u5 + e ) 2 ]} 



— (2i arctan[2(w - e )] + 2z arctan[2(w + e )] - ln[l + 4(S - £ ) 2 ] - ln[l + 4(5 + e ) 2 ] + 2 ln[l + 4eg])l , (B8) 



and 



n<(S5) 



X 



-{2a; arctan[2(a; — e )] + 2cj arctan[2(w + e )] 



2'K LO + U> 3 

+ ln[l + 4(2 - e ) 2 ] + ln[l + 4(2 + e ) 2 ] - 2 ln[l + 4^]}fl(-3). 
Clearly, ImII r (a;) and II < (w) yield a finite lifetime for the vibrational mode. 



(B9) 



(g) 



o 



FIG. 6: Top: The thick dotted line denotes the Green function of the 
vibrational mode with the RPA term, see Eqs. dB It and dB2t . Bot- 
tom: The diagrams required when the RPA term is included besides 
the diagrams in Fig. Q][with the replacement of the dotted line in Fig. 
[Tjby the thick dotted line, as shown in the top part]. 



* Electronic address: akiko@bgu.ac.il 

* Also at Tel Aviv University, Tel Aviv 69978, Israel 

1 M. A. Reed, C. Zhou, C. J. Muller, T. P. Burgin, and J. M. Tour, 
Science 278, 252 (1997); J. Chen, M. A. Reed, A. M. Rawlett, 
and J. M. Tour, ibid. 286, 1550 (1999). 

2 J. Reichert, R. Ochs, D. Beckmann, H. B. Weber, M. Mayor, and 
H. v. Lohneysen, Phys. Rev. Lett. 88, 176804 (2002). 

3 N. B. Zhitenev, H. Meng, and Z. Bao, Phys. Rev. Lett. 88, 226801 
(2002). 

4 S. Kubatkin, A. Danilov, M. Hjort, J. Cornil, J. Bredas, N. Stuhr- 
Hansen, P. Hedegard, and T. Bjornholm, Nature (London) 425, 
698 (2003). 

5 J. G. Kushmerick, J. Lazorcik, C. H. Patterson, R. Shashidhar, D. 
S. Seferos, and G. C. Bazan, Nano Lett. 4, 639 (2004). 

6 X. H. Qiu, G. V. Nazin, and W. Ho, Phys. Rev. Lett. 92, 206102 
(2004). 

7 H. Park, J. Park, A. K. L. Lim, E. H. Anderson, A. P. Alivisatos, 
and P. L. McEuen, Nature (London) 407, 57 (2000); A. N. Pa- 



supathy, J. Park, C. Chang, A. V. Soldatov, S. Lebedkin, R. C. 
Bialczak, J. E. Grose, L. A. K. Donev, J. P. Sethna, D. C. Ralph, 
and P. L. McEuen, Nano Lett. 5, 203 (2005). 
B. J. LeRoy, S. G. Lemay, J. Kong, and C. Dekker, Nature (Lon- 
don) 432, 371 (2004). 

R. H. M. Smit, Y. Noat, C. Untiedt, N. D. Lang, M. C. van Hemert, 
and J. M. van Ruitenbeek, Nature (London) 419, 906 (2002); D. 
Djukic, K. S. Thygesen, C. Untiedt, R. H. M. Smit, K.W. Jacob- 
sen, and J. M. van Ruitenbeek, Phys. Rev. B. 71, 161402 (2005); 
W. H. A. Thijssen, D. Djukic, A. F. Otte, R. H. Bremmer, and J. 
M. van Ruitenbeek, Phys. Rev. Lett. 97, 226806 (2006). 
S. Sapmaz, P. Jarillo-Herrero, Ya. M. Blanter, C. Dekker, and H. 
S. J. van der Zant, Phys. Rev. Lett. 96, 026801 (2006). 
O. Tal, M. Krieger, B. Leerink, and J. M. van Ruitenbeek, Phys. 
Rev. Lett. 100, 196804 (2008); M. Kiguchi, O. Tal, S.Wohlthat, F. 
Pauly, M. Krieger, D. Djukic, J. C. Cuevas, and J. M. van Ruiten- 
beek, Phys. Rev. Lett. 101, 046801 (2008). 
R. Leturcq, C. Stampfer, K. Interbitzin, L. Durrer, C. Heirold, E. 



10 



Mariani, M. G. Schultz, F. von Oppen, and K. Ensslin, Nature 
Physics 5, 327 (2009). 

13 M. Galperin, M. A. Ratner, and A. Nitzan, Nano Lett. 4, 1605 
(2004). 

14 A. Mitra, I. Aleiner, and A. J. Millis, Phys. Rev. B 69, 245302 
(2004). 

15 D. A. Ryndyk and J. Keller, Phys. Rev. B 71, 073305 (2005). 

16 O. Hod, R. Baer, and E. Rabani, Phys. Rev. Lett. 97, 266803 
(2006); J. Phys: Condens. Matter, 20, 383201 (2008). 

17 R. Egger and A. O. Gogolin, Phys. Rev. B 77, 1 13405 (2008). 

18 R-P. Riwar and T. L. Schmidt, Phys. Rev. B 80, 125109 (2009). 

19 O. Entin-Wohlman, Y. Imry, and A. Aharony, Phys. Rev. B 80, 
035417 (2009); O. Entin-Wohlman, Y. Imry, and A. Aharony, 
Phys. Rev. B 81, 1 13408 (2010). 

20 J. Koch and F. von Oppen, Phys. Rev. Lett. 94, 206804 (2005); J. 
Koch, M. Semmelhack, F. von Oppen, and A. Nitzan, Phys. Rev. 
B 73, 155306 (2006). 

21 L. Muhlbacher and E. Rabani, Phys. Rev. Lett. 100, 176403 
(2008). 

22 M. Galperin, A. Nitzan, and M. A. Ratner, Phys. Rev. B 74, 
075326 (2006). 



23 T. L. Schmidt and A. Komnik, Phys. Rev. B 80, 041307(R) (2009). 

24 R. Avriller and A. Levy Yeyati, Phys. Rev. B 80, 041309(R) 
(2009). 

25 F. Haupt, T. Novotny, and W. Belzig, Phys. Rev. Lett. 103, 136601 
(2009); Phys. Rev. B 82, 165441 (2010). 

26 D. F. Urban, R. Avriller, and A. Levy Yeyati, Phys. Rev. B 82, 
121414 (2010). 

27 R. Avriller, lai-Xiv: 1007.44501 

28 J.-X. Zhu and A. V. Balatsky, Phys. Rev. B 67, 165326 (2003). 

29 S. Maier, T. L. Schmidt, and A. Komnik. larXiv: 10 10.29 1 8 1 

30 B. Kubala and F. Marquardt, Phys. Rev. B 81, 1 15319 (2010). 

31 L. V. Keldysh, Zh. Ekps. Teor. Fiz. 47, 1515 (1964) [Sov. Phys. 
JETP 20,1018 (1965)]. 

32 H. J. W. Haug and A-P. Jauho, Quantum Kinetics in Transport and 
Optics of Semiconductors (Springer Verlag 2008). 

33 A. Ueda, O. Entin-Wohlman and A. Aharony, Nanoscale Res. 
Lett, (in press). 

34 G. D. Mahan, Many-Particle Physics (Plenum Press, New York, 
1990). 



